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Magnetic pyrochlore oxides, including the spin ice materials, have proved to be a rich field for the study of geomet¬ 
rical frustration in 3 dimensions. Recently, a new family of magnetic oxides has been synthesised in which the half of 
the tetrahedra in the pyrochlore lattice are inflated relative to the other half, making an alternating array of small and 
large tetrahedra. These “breathing pyrochlore” materials such as LiGaCr408, LiInCr408 and Bas Yb2Zn50ii provide new 
opportunities in the study of frustrated magnetism. Here we provide an analytic theory for the ground state phase dia¬ 
gram and spin correlations for the minimal model of magnetism in breathing pyrochlores: a classical nearest neighbour 
Heisenberg model with different exchange coefficients for the two species of tetrahedra. We find that the phase diagram 
comprises a Coulombic spin liquid phase, a conventional ferromagnetic phase and an unusual antiferromagnetic phase 
with lines of soft modes in reciprocal space, stabilised by an order-by-disorder mechanism. We obtain a theory of the 
spin correlations in this model using the Self Consistent Gaussian Approximation (SCGA) which enables us to discuss 
the development of correlations in breathing pyrochlores as a function of temperature, and we quantitatively characterise 
the thermal crossover from the limit of isolated tetrahedra to the strongly correlated limit of the problem. We compare 
the results of our analysis with the results of recent neutron scattering experiments on LiInCr408. 


1. Introduction 

Systems exhibiting geometrical frustration have proved 
over recent years to be fertile ground for the discovery 
of novel emergent phenomena.Amongst frustrated mag¬ 
netic systems there are now several known materials whose 
low temperature physics displays novel collective behaviour 
which cannot be understood using the usual Landau picture, 
based on broken symmetries. 

At the forefront of the research interest in these systems 
has been the study of magnetism on the pyrochlore lattice^ -a 
network of corner sharing tetrahedra. Nearest neighbour an¬ 
tiferromagnetic coupling between spins on this lattice fails to 
select a unique ordered state and leads to an exponentially 
large ground state manifold. This results in the forma¬ 
tion of a classical spin liquid with extensive residual entropy 
and algebraically decaying spin correlations.^^’The ground 
state manifold of this spin liquid is defined by a flux conser¬ 
vation law V • B = 0 and the excitations which violate this 
law are deconflned “charges” which interact via an emergent 
Coulomb law. For this reason, this spin liquid is known as a 
“Coulombic spin liquid”. 

The classical Coulomb spin liquid flnds realization in sev¬ 
eral experimental systems. The spin ice materials R2M2O7 
(R=Ho, Dy; M=Ti, Sn) at low temperature are famous ex¬ 
amples of Coulomb liquids of Ising spins.Meanwhile, the 
chromium spinels MCr204 (M=Zn, Cd, Hg) can be under¬ 
stood in terms of a Coulomb liquid of vector spins interacting 
with fluctuations of the lattice. As a consequence these mate¬ 
rials have emerged as canonical examples of the phenomenon 
of “order-by-distortion” in which magnetic frustration is re¬ 
lieved by a structural transition.Calculations based on a 
classical effective model for the spins with coupling to the 
lattice predict a rich phase diagram in an applied magnetic 
fleld^^“^^ which was observed in a series of remarkable exper¬ 
iments at fields up to 600 T.^^“^^ 



Fig. 1. (Color online) The breathing pyrochlore lattice. Like the pyrochlore 
lattice, the breathing pyrochlore lattice is composed of corner-sharing tetra¬ 
hedra. However, in the breathing pyrochlore lattice these tetrahedra alternate 
in size. The spins reside at the vertices of this lattice. The two species of tetra¬ 
hedra (here labelled ‘A’ and ‘B’, colored black and red) have different bond 
lengths Ja and ds and different exchange coefficients Ja and Jb. 


Recently, a new class of magnetic materials has been syn¬ 
thesised which promises an interesting variation on this pic¬ 
ture. In the spinels LiInCr408 and LiGaCr408^^“^^ the 

Cr^+ ions, carrying spin S = reside on a “breathing” py¬ 
rochlore lattice depicted in Fig. 1 . In this lattice the tetrahedra 
of the pyrochlore lattice are alternately small and large, with 
the consequence that the two sets of tetrahedra have differ¬ 
ent exchange coefficients. Similarly, in the effective S = \ 
magnet, Ba3Yb2Zn50ii, the magnetic Yb^+ ions reside on a 
breathing pyrochlore lattice.^'^ 

As with Cr spinels, these oxide materials have antiferro¬ 
magnetic nearest-neighbour exchange interactions.^^ How¬ 
ever this interaction has contributions from both direct and 
indirect exchange, with opposite sign, and is very sensitive 
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small changes in bond angles or distances^^’^^ . As a re¬ 
sult, there is a very real possibility that breathing pyrochlores 
based on larger anions, such as S or Se, might have ferro¬ 
magnetic interactons, in at least one of the the two sets of 
tetrahedra. 

Theoretical studies of the pyrochlore magnets with un¬ 
equal exchange interactions on the two species of tetrahe¬ 
dra have concentrated on antiferromagnetic interactions, in 
the quantum limit S = 1 / 2 , as a route to understanding the 
S = 1/2 Heisenberg model on a pyrochlore lattice.Per¬ 
forming perturbation theory about the limit of independent 
independent tetrahedra, Harris et al?^ and Tsunetsugu"^^ have 
explored the possibility of a novel dimerized ground state. 
Meanwhile Canals and Lacroix,^^’^^ considered the possibil¬ 
ity of a quantum spin liquid. At present, there is no established 
study of the Heisenberg model on the breathing pyrochlore 
lattice in the presence of ferromagnetic interactions, or treat¬ 
ment of spin correlations within the high-temperature param¬ 
agnetic phase. The recent synthesis of breathing pyrochlore 
magnets makes this a pressing issue. 

In this article we provide a comprehensive treatment of the 
ground state phase diagram and spin correlations of the near¬ 
est neighbour Heisenberg model on the breathing pyochlore 
lattice, for all signs and ratios of exchange couplings. The 
classical treatment should be particularly relevant for the rel¬ 
atively large-spin 5 = 1, Cr^+ systems, in the temperature 
regime above the structural transitions which are observed in 
LiInCr408 at 15.9 K and LiGaCr408 at 13.8 

We find that the nature of the classical ground state de¬ 
pends only the signs of the two exchange coefficients asso¬ 
ciated with the two species of tetrahedra, not on their relative 
magitudes. Where both exchange coefficients are antiferro¬ 
magnetic the system enters a Coulomb phase at low tempera¬ 
ture. Where both exchange coefficients are ferromagnetic the 
system orders into a simple ferromagnetic phase. Where the 
exchange coefficients differ in sign, the classical ground state 
manifold has a degeneracy of C)(L), where L is the linear di¬ 
mension of the system, with antiferromagnetic [001] planes 
becoming effectively decoupled from one another. This deg- 
neracy is then susceptible to being lifted by fluctuations, sta¬ 
bilising antiferromagnetic order, meaning that this region of 
the phase diagram combines the properties of dimensional re¬ 
duction and order-by-disorder noted in other pyrochlore mag- 
nets.^^"^^ By using the Self Consistent Caussian Approxima¬ 
tion (SCCA)'^^“'^^ we are then able to calculate the behaviour 
of the spin correlations as a function of the exchange parame¬ 
ters and temperature. 

The article is structured as follows: 

In Section 2 we obtain the phase diagram of the classical 
Heisenberg model on the breathing pyrochlore lattice. 

In Section 3 we discuss the case of ferromagnetic coupling. 
We use the SCCA to show the temperature development of the 
spin correlations and quantify the thermal crossover from iso¬ 
lated tetrahedra to long range correlations by deriving equa¬ 
tions for the temperature dependence of the ferromagnetic 
correlation length. 

In Section 4 we discuss the case of exchange parameters 
with differing signs. We show that the degenerate zero-energy 
modes which occur in this case lead to a “square ring” struc¬ 
ture observable in the neutron scattering response. 

In Section 5 we discuss the case of antiferromagnetic cou- 


Jb 

J 



Fig. 2. (Color online) Classical ground state phase diagram of the breath¬ 
ing pyrochlore lattice Heisenberg model [Eq. (1)]. The ground state mani¬ 
fold is determined by the signs of the two couplings Ja and Jb- Where both 
couplings are ferromagnetic than the ground state is a simple collinear fer- 
romagnet. Where both couplings are antiferromagnetic the system enters a 
macroscopically degnerate Coulomb phase, in which the ground state mani¬ 
fold contains all states where every tetrahedron has vanishing magnetisation. 
If Ja and Jb have opposite signs than the set of classical ground states is 
described by decoupled [001] planes and has a degeneracy 0{L). This de¬ 
generacy will be lifted by the order-by-disorder mechanism. Due to a cor¬ 
respondence with Heisenberg model on the FCC lattice [Eq. (11)] we expect 
this to favour antiferromagnetic states with ordering vector qord = (1,0,0) 
[Fig. 3]. The green points indicate the parameter sets for which we plot the 
spin correlation functions in Figs. 4-9. 


pling which leads to a Coulombic spin liquid at low tempera¬ 
ture. We show how the algebraic correlations of the spin liquid 
emerge out of the short range, single tetrahedron, correlations 
with decreasing temperature. 

In Section 6 we conclude with a summary of our re¬ 
sults and a comparison of our predictions with the results 
of neutron scattering experiments on the breathng pyrochlore 
LiInCr408.^^ 

2. Ground state phase diagram 

We consider the Heisenberg Hamiltonian on the breathing 
pyrochlore lattice 

-Hbreathing = -/a ^ S,- • S,- + 7^2 S,- • S,- ( 1 ) 

<0>A <0>B 

The two terms in Eq. ( 1 ) correspond to the two sets of tetra¬ 
hedra on the pyrochlore lattice, which we label A and B 
[cf. Fig. 1 ]. Every bond of the lattice belongs uniquely to one 
tetrahedron and is associated with an exchange coupling Ja 
or Jb depending on which set of tetrahedra it belongs to, as 
illustrated in Fig. 1 . 

The classical ground state phase diagram of Eq. ( 1 ) is 
shown in Fig. 2 . As is commonplace for two parameter mod¬ 
els, we can re-express the exchange coefficients in terms of an 
overall energy scale / and an angle 

Ja = /cos(?^), Jb = /sin(?^) (2) 

There are three distinct phases on the phase diagram: 


2 





J. Phys. Soc. Jpn. 


FULL PAPERS 


(1) For 0 < ?^ < I (/a > 0,> 0) the ground state 
manifold contains all states where the magnetisation of 
every tetrahedron in the lattice vanishes and the system 
is in a Coulomb phase at low temperature. 

(2) For —TT < ?^ < — I (/a < 0, < 0) the ground state is 

a collinear ferromagnetic state 

(3) For — I < 1 ^ < 0 (Ja > 0,Jb < 0) or | < § < n 
(Ja < 0,Jb > 0) the ground state is comprises an 0{L) 
set of antiferromagnetic states with wavevectors of the 
form of q <5 = (d, 0,1) (and those related by symmetry). 
This ground state manifold may be mapped directly to 
the classical ground state manifold of the nearest neigh¬ 
bour FCC antiferromagnet,^^’^^ as we shall show below. 

The phase diagram shown in Fig. 2 may be derived us¬ 
ing the following argument. As in the case of the Heisenberg 
model on the ordinary (non-breathing) pyrochlore lattice, the 
Hamiltonian [Eq. (1)] can be rewritten as a function only of 
the magnetisation of each tetrahedron. We define, for a tetra¬ 
hedron t, the magnetisation density 

(3) 

iEt 

where the sum runs over the four spins at the vertices of tetra¬ 
hedron t. Eq. (1) can then be rewritten solely in terms of the 
variables 

-Hbreathing = 

teA teB 


{Ja + Jb) 

2 


NS^ 


( 4 ) 


It is clear from Eq. (4) that when Ja and Jb are both anti¬ 
ferromagnetic, any state where 


= 



Mt 


( 5 ) 


is a classical ground state. This constraint gives rise to a 
Coulomb phase. 

It is similarly clear from Eq. (4) that when Ja and Jb are 
both ferromagnetic, the ground state is obtained when every 
tetrahedron has a maximal value of 

= 1 V r ( 6 ) 


This condition is ony satisfied by collinear ferromagnetic 
states. 

Where the signs of Ja and Jb differ we can see from Eq. (4) 
that the condition to be in a ground state is that all ferromag¬ 
netic tetrahedra have aligned spins but the intervening antifer¬ 
romagnetic tetrahedra have vanishing magnetisation. Taking, 
without loss of generality 


/a > O? Jb ^ ^ 



(a) The FCC lattice as a network of edge¬ 
sharing tetrahedra 



(b) The (100) state on the FCC lattice 



(c) The (100) state on the breathing py¬ 
rochlore lattice 

Fig. 3. (Color online). The mapping to the FCC antiferromagnet when Ja 
and Jb have opposite signs and the resulting antiferromagnetic order. Each 
sublattice of tetrahedra (A and B) of the breathing pyrochlore lattice forms 
an FCC lattice, which can be viewed as a network of edge sharing tetrahe¬ 
dra, as seen in A. Where the exchange coefficients Ja and Jb differ in sign, the 
classical ground state manifold of Eq. (1) is the same as the ground state man¬ 
ifold of the nearest neighbour antiferromagnet on the FCC lattice, with each 
ferromagnetic tetrahedron of the breathing pyrochlore model corresponding 
to a single spin of the FCC model. This model has an 0{L) set of ground 
states with fluctuations favouring collinear states with ordering wavevector 
qord = (1,0,0) (and states related by global symmetries). As a result, the 
FCC model orders into the configuration shown in (b). Since the two models 
have an exactly corresponding structure of soft modes, the order by disorder 
mechanism should favour the analagous state on the breathing pyrochlore 
lattice [shown in (c)] when Ja and Jb have opposite signs. 


we require 

= OV f G A (7) 

= IV f G 5 ( 8 ) 

These constraints lead to an equivalent classical ground 
state manifold to the nearest neighbour Heisenberg model on 
the FCC lattice. We can see this by considering the Hamilto¬ 


nian [Eq. (1)] restricted to the subspace where Eq. ( 8 ) is per¬ 
fectly obeyed. Since this fixes the exchange energy of all the 
B tetrahedra, the Hamiltonian within this set of configurations 
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becomes 

(9) 

teA 

Since all the spins on each B tetrahedron are ferromagneti- 
cally aligned the magnetisation of the A tetrahedra is deter¬ 
mined uniquely by the total magnetisation of the four sur¬ 
rounding B tetrahedra 

= \ Yj 

tB^n.n.tA 

Inserting Eq. (10) into Eq. (9) and using Eq. (8) we obtain 
a nearest neighbour antiferromagnetic interaction between A 
tetrahedra. 

n™a^hing = JaS^ m,, • m,, + (11) 

where the sum nearest neighbours of the 

ECC lattice of B tetrafiedra. 

The ECC lattice may be depicted as a network of edge¬ 
sharing tetrahedra as shown in Eig. 3(a). A ground state for 
the nearest neighbour antiferromagnet on this lattice is ob¬ 
tained when each of these edge-sharing tetrahedra have van- 
shing magnetisation. Under the mapping from the breathing 
pyrochlore model [Eq. (4)] to the ECC model [Eq. (11)] this 
requirement is equivalent to Eq. (7). 

One ground state configuration of Eq. (11) is depicted in 
Eig. 3(b). This is an antiferromagnetic state with ordering 
wavevector qord = (1,0,0). The corresponding state on the 
breathing pyrochlore lattice is depicted in Eig. 3(c) and one 
can readily verify that it satisfies Eqs.(7)-(8). 

Eurther ground states can be generated from the configu¬ 
ration in Eig. 3(b), by taking (e.g.) [001] planes of spins and 
rotating them with respect to the rest of the lattice. This oper¬ 
ation keeps the system in the classical ground state manifold. 
By generating ground states using these operations we see that 
the ground state manifold contains states with wavevectors 

qs = {h0,S) (12) 

(and those wavevectors related to Eq. (12) by lattice symme¬ 
tries). 

Erom these considerations we can see that states with 
wavevectors 

qord = (1,0,0) (13) 

are special because they live at the junctions of the ground 
state manifold, where two lines of classical ground states 
(in the case of Eq. (13) those with qs = (1,0, d) and 
qs = (l,d, 0)) meet. The resulting freedom to fluctuate for 
zero energy cost should favour these states for selection via 
the order-by-disorder mechanism. 

In the nearest neighbour ECC model it is known that 
fluctuations do indeed favour states with ordering wavevec¬ 
tor qord = (1,0,0) (referred to in the literature as Type 
I AEM states).This entropic preference for states with 
Qord = (1,0,0) results in an first order phase transition into 
these states.^^ 

Since the ground state manifold of the breathing pyrochlore 
model [Eq. (1)] with Ja < 0, Jb > 0 (or equivalently Ja > 
0 , Jb < 0) is the same as in the ECC model [Eq. (11)], the 


structure of the soft modes is also the same, and it therefore 
seems reasonable to postulate that the ground state selection 
will also be the same. 

In the limit \Ja\ » \Jb\, for J a <0, Jb >0, i.e. where the 
ferromagnetic exchange is much stronger than the antiferro¬ 
magnetic, fluctuations which are internal to the ferromagnetic 
tetrahedra will cost a lot of energy and [Eq. (11)] will be a 
faithful treatment not only of the ground state manifold but 
also of the low energy fluctuations. Therefore, in this limit at 
least, the ground state selection must be the same as in the 
ECC model. Since tuning away from the limit \Ja\ » \Jb\ 
does not change the structure of the ground state manifold we 
make the conjecture that the qord = (1,0,0) state is selected 
across the whole region where Ja and Jb have opposite signs, 
although this would need to be confirmed by simulation. 

In the temperature range above the phase transition, two 
planes neighbouring [001] planes with short range order at 
q = (1,0,0) will become effectively decoupled and fluctu¬ 
ate essentially independently. This gives rise to an apparent 
dimensional reduction, which will manifest itself as line like 
features in the structure factor in the paramagnetic phase, in 
a similar phenomenon to that observed in the paramagnetic 
phase of the rare earth pyrochlore Yb 2 Ti 207 .^^’'^^’'^^ We will 
discuss this point further where we calculate the spin correla¬ 
tions in this region of the phase diagram in Section 4. 

3. Ja < 0, Jb < 0: ferromagnetic phase 

In this Section we analyse the spin correlations where the 
exchange coefficients Ja and Jg are ferromagnetic. Here, and 
in the following two Sections where we analyse the other re¬ 
gions of the phase diagram, we use the Self Consistent Gaus¬ 
sian Approximation (SCGA) to calculate the spin correlation 
functions. This approach is known to be qualitatively succes- 
ful in treating the spin correlations of classical vector spin 
models on the pyrochlore lattice. 

In this treatment of the spin correlations, the “hard spin” 
constraint that each vector spin have fixed length 

S,-S,=5^ V i (14) 

is relaxed such that Eq. (14) is satisfied only on average: 

(SrSi} = S^ V i (15) 

Eq. (15) is enforced by means of a Lagrange multiplier A. 
Specifically, we write 

^"Hbreathing ^ ^"^breathing 

= y8<Hbreathing + (4 - ySeg) XI (16) 

i 

Here p is the inverse temperature and 


where Eq is the ground state energy. In Eq. (16) we have sub¬ 
tracted the constant term eq from the Hamiltonian so 
is positive definite as long as T > 0. This step 
makes no difference to the results of our analysis, but is math¬ 
ematically convenient. 

We define Eourier transformed variables 
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Fig. 4. (Color online) Spin structure factor <S(q) [Eq. (??)] in the Ferromagnetic region of the phase diagram of TTbreathing [Eq. (1)] for (close 

to the non-breathing limit Ja = Jb) and (close to the isolated tetrahedron limit Ja = 0) [cf. Fig. 2]. <S(q) is calculated from the Self Consistent 

Gaussian Approximation (SCGA) described at the beginning of Section 3. With decreasing temperature the scattering sharpens around q = (0,0,0) and 
other ferromagnetic zone centers. The sharpening is less pronounced at equal temperature for § = because the growth of the correlation length at low 
temperatures is controlled by the product JaJb which vanishes in the isolated tetrahedron limit [Eqs. (30)-(36)]. 


where a = x,y,z indexes the spin components, i = 0,1,2, 3 
indexes the four sublattices of the breathing pyrochlore lattice 
and Yiv belonging to sublattice i. 

Re-writing yST^breShing i" terms of 5“(q) [Eq. 

(18)] we obtain 


^breathing = ^ 2 2 2 5“(-q)Mi(q)5“(q). 

q u ij 

The components of 4 x 4 matrix Af; 7 (q) are 

Mijiq) 


(19) 


— 26ij(A — pet)) -|- 


Ja exp ( (e,- - ey) 


where 


+Jb exp 



/O 1 1 1\ 

110 1 
(l 1 1 o) 

eo = (1,1,1) ei = (1,-1,-1) 
62 = (-1,1,-!) 63 = (-1,-1,!) 


( 20 ) 


( 21 ) 


( 22 ) 


and Ja, Js are the bond lengths on the A and B tetrahedra 
respectively [Fig. 1]. 
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(a) ^ = 


-Un . r 
16 ’ / 


= 2(a) 



The partition function within the SCGA is then 

zscGA=J n n (-^t 


SCGA \ 
breathing j 


(23) 


Since is Gaussian the spin correlation functions 

within the SCGA are given by 


<5r(-q)5Kq)> = <5'^’'[M(q)-‘],r 


(24) 


The Lagrange multiplier A is then determined self consistently 
from Eqs. (15) and (24). 

The results of this analysis, for ferromagnetic exchange pa¬ 
rameters Ja,Jb<^ are shown in 4 and Fig. 5. 

Fig. 4 shows the spin structure factor 


*S(q)=2<S,-{-q)-SXq)> 


(25) 


IJ 


in the plane q = {h, h, 1) . The structure factor is shown for 
two temperatures 7=5 and 7 = 2 for ^ (close to 

the “non-breathing” limit J a = Jb) and 6 = ^ (close to the 
isolated tetrahedra limit /a = 0 )- Fig- 5 shows the angle inte¬ 


gral of Eq. (25), for the same sets of exchange parameters and 
a range of temperatures, as would be measured in a neutron 
scattering experiment on a powder sample 


<^pow(2) — 


^ r d6d(p sin{6)S (q) (26) 

An J 


q = 2(sin(^) cos( 0 ), sin(^) sin( 0 ), cos{6)) (27) 

In both cases there is a sharpening of the scattering around 
into peaks around q = (0,0,0) (and other Brillouin zone cen¬ 
ters associated with ferromagnetic Bragg peaks) as the tem¬ 
perature approaches the ordering temperature Tfm- Within 
the SCGA this transition occurs when A{T) = 0. Since the 
SCGA includes unphysical fluctuations which modulate the 
spin length, the resulting estimate of the transition temper¬ 
ature is likely to be an underestimate of the true transition 
temperature. 

The sharpening of the scattering is less pronounced (for 
equal temperature) where the exchange parameters are closer 
to the isolated tetrahedra limit Ja = 0. In essence, where there 
is a large disparity between Ja and Jb the correlation length 
is limited by the weaker of the two exchange parameters. 

This can be understood quantitatively by considering the 
behaviour of the ferromagnetic correlation length, obtained 
by expanding the SCGA calculation of Eq. (25) around 
q = (0,0,0). For this purpose, and for the calculations in Sec¬ 
tions 4 and 5, we will set the ratio of bond lengths on the 
tetrahedra [Fig. 1] equal to unity 


2n 

(b) 0 = ^; 1= 5(»), I = 4(,), I = 3(»), I = 2(A) 

Fig. 5. (Color online) Angle integrated structure factor <Spow(G) [Eq. 
(27)] in the ferromagnetic region of the phase diagram [Fig. 2] of 
Tfbreathing [Eq. (1)]. The Structure factor is calculated from the SCGA as de¬ 
scribed at the beginning of Section 3. Calculation is shown for and 

. Where the disparity between Ja and Jb is greater the feromagnetic 
correlations build up more slowly with decreasing temperature. 


^ — I 

dB 


(28) 


This assumption makes the analysis more transparent, and 
is not unphysical since in the real breathing pyrochlore ma 
terials ^ is very close to unity, e.g. 0.95 in LiInCr 408 . 
Where we come to compare our calculations to experiments 
on LilnCraOs in Section 6 we restore the true value of — 


30 


r4^« AXX V. VV^ X^O.V^X^ .XX^ .X«^ VCXXUA^ V.X . 

In the region approaching q = (0,0,0) we have 


'S(q) = 


1 


KT) (1+^? 


+ 0{q% 


(29) 


eff, 


The square of the ferromagnetic correlation length may 
be written as the sum of two terms, with different behaviours 
approaching the ordering transition 


^eff ~ ^0 + ^FM* 


These two terms are 


= (t)’ 
^ (?) 


P\Ja + Jb\ 


2{A{T)+2I3\Ja + Jb\) 

2 A0^JaJbA{T) 


(30) 


(31) 


(32) 


A{TY{A{T) + 2I3\Ja + Jb\) 

where gq is the cubic lattice parameter of the breathing py 
rochlore lattice. 

Approaching the ordering transition at T = Tfm, we have 


lim A(T) = 0. 


(33) 


In this limit the correlation length saturates to a length scale 
of the order of the nearest neighbour distance 




(34) 
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Fig. 6. (Color online) Spin structure factor <S(q) [Eq. (25)] for TYbreathing [Eq. (1)] in the region of the phase diagram [Fig. 2] where Ja and Jb have opposite 
signs, in the plane q = (h, k, 0). For all values of ^ in this region the correlations evolve towards a “square ring” pattern, which comes from interesecting lines 
of soft modes with wavevectors of the form q = (1, d, 0) 
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Meanwhile the correlation length diverges as j 

AP^JaJbKT) 1 


lim^FMW^ = (^) 


Therefore as T 


A{TY2P\Ja + Jb\ 

T FM we have 

^eff ^FM 

In the high temperature limit yS ^ 0 we have 

_ P\Ja + Jb\ 

^ U1 Cn — I 


^FM » ^0 


= (?) 
= (f) 


A{T) 

and so at high temperatures we have 


^P^JaJb _ ^2 




FM 




(35) 

(36) 

(37) 

(38) 

(39) 


Seen in this light, the two correlation lengths and ^fm are 
a quantitative measure of the extent to which the system may 
be thought of as a set of isolated tetrahedra. The crossover 
from the high temperature limit described by Eq. (39) to the 
strongly correlated regime described by Eq. (36) is a ther¬ 
mal crossover from the physics of isolated tetrahedra to the 
physics of a long-range correlated ferromagnet. The relative 
magnitudes of and ^fm give a criterion for deciding whether 
a given ferromagnetic breathing pyrochlore at a given temper¬ 
ature is in the “isolated tetrahedron” limit. 

When either Ja or Jp vanishes, ^fm also vanishes for all 
temperatures, since the system is then exactly described by a 
set of isolated tetrahedra. 


4. J A and Jb with opposite signs: dimensional reduction 
and order by disorder 

We now turn to discuss the case where Ja and Jb have op¬ 
posite signs. Eor concreteness we will take throughout this 
section 


Ja >0, Jb< 0. (40) 

As discussed in Section 2, the ground state manifold in 
this case has a degeneracy of C)(L), with states in the ground 
state manifold having associated wavevectors of the form 
= (1,0,^). At some finite temperature Tobd the system 
will order via the order-by-disorder mechanism. The rela¬ 
tionship between this model and the classical antiferromagnet 
on the ECC lattice lEq. (11)] suggests that this will occur via 
a first-order phase transition^ ^ and that the resulting ordered 
state will have wavevector qord = (1,0,0) (or those related by 
symmetry). 

The development of these correlations in the plane 
q = (h, k, 0) as a function of temperature and lEq. (2)] is 
shown in Eig. 6, for “f ’ cases, the 

scattering develops with decreasing temperature towards the 
“square ring” pattern visible in Eigs. 6 (g)-(i), which is formed 
from intersecting lines of scattering associated with the soft 
modes at qs = (1, 0), modulated by a zone dependent form 

factor. These bright lines of scattering signal an effective de¬ 
coupling of neighbouring [001] planes- which is somewhat 
reminiscent of the dimensional reduction observed in the rare- 
earth pyrochlore Yb2Ti207.'^^’'^^’'^^ In that system neighbour¬ 
ing [111] planes become effectively decoupled and the conse¬ 
quence is rods of scattering along the {h, h, h) directions of re- 



(a)^ = ^ 5 - 5(#), t = 2(-),5 = 1W,5 = 
0.5(a), t = o.2( ) 



(b) = ^; £ = 5(»), £ = 2{m)A = !(♦), 5 = 

0.5(A), Z = 0.2( ) 



In 

(c)i>= Zj = 5(«), z = 1(A) Z = 1(,),F ^ 

0.5(a), Z = o.2( ) 

Fig. 7. (Color online) Angle integrated structure factor <Spow (G) [Eq. (27)] 
for T/breathing [Eq. (1)] in the region of the phase diagram [Fig. 2] where Ja 
and Jb have opposite signs. At high temperatures, correlations are limited 
to single tetrahedra and are ferromagnetic or antiferromagnetic depending 
on the sign of the dominant exchange paramter. At low temperatures these 
correlations evolve into an increasingly sharp peak near ^ = 1, which is 
the result of antiferromagnetic correlations between internally ferromagnetic 
tetrahedra. The peak near ^ = 1 comes from integrating the square ring 
structure seen in Fig. 6 (g)-(i). 


ciprocal space. Here, the rods are along the (1,0, /z) directions 
and are modulated by a zone dependent form factor resulting 
in the square ring structure. 

In the case ^ [Eigs. 6 (b), (e), (h)], where \Ja\ = \Jb\ 
the scattering is only very weakly q dependent at j =5. This 
reflects an effect by which the frustration arising from the op¬ 
posite exchange parameters, effectively reduces the temper¬ 
ature at which appreciable q dependence in the correlations 
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arises. The “square ring” structure arises smoothly out of this, 
becoming increasingly sharp as the temperature is lowered. 

In the case ^ [Figs. 6 (a), (d), (g)], the antiferromag¬ 
netic Ja has a much larger magnitude than the ferromagnetic 
/fi. The high temperature correlations therefore simply re¬ 
flect antiferromagnetic nearest-neighbour correlations on the 
A tetrahedra, with a broad minimum around q = (0,0,0) and 
a broad minimum at the neighbouring Brillouin zone centers 
(e.g.)q = (2,0,0). 

At a flnite temperature given by the equation 

t^AFM i^T^ring ^ — 6/5(6 /a “ 2/^) 

^ring ^(^AFm) = -- (41) 

the maximum at zone centers neighbouring the first Bril¬ 
louin zone (e.g. q = (2,0,0)) inverts into a minimum, and 
a set of continuous set of maxima appear appear between 
q = (0,0,0) and the neighbouring Brillouin zone centers. 
This set of maxima evolves in to the “square ring” feature 
seen at low temperature. This temperature evolution of the 
correlations can be seen in angle integrated scattering shown 
in Fig. 7(a), where at broad maximum centred at 



at high temperature turns into a minimum, with an increas¬ 
ingly sharp peak evolving towards 

Qcio _ ^ 

In 

The peak in the powder integral always retains a flnite width 
and is shifted slightly to the right of ^ = 1 because it has a 
contribution from all wavevectors of the form q = (1,0, d). 

A similar picture describes the scattering when the ferro¬ 
magnetic exchange parameter dominates the antiferromag¬ 
netic one. This is the case when § = ^ [Figs. 6 (c), (f), 
(i)]. Here the high temperature data has a broad maximum 
at q = (0,0,0) reflecting ferromagnetic nearest neighbour 
correlations on the B tetrahedra. This maximum inverts at a 
temperature given by 

t^FM — '^Ja) 

(42) 

and a continuous set of maxima appear which evolve towards 
the “square ring” feature seen at low temperature. The same 
evolution can be seen in the powder integrated scattering in 
Fig. 7(c). 

It is interesting to note that a “square ring” feature, rather 
similar to that seen in Fig. 6 has also been seen in reverse 
Monte Carlo reconstructions of the scattering from the dou¬ 
ble perovskite material Ba 2 YRu 06 in which the Ru^+ ions 
form an FCC lattice.This material orders in to q = (1,0,0) 
antiferromagnetic state shown in Fig. 3(b) below T = 36K. 

5. Ja > 0, Jb > 0: Coulombic spin liquid 

We now turn to discuss the spin correlations in the anti¬ 
ferromagnetic region of the phase diagram where the system 
enters a Coulomb phase at low temperature. These are shown 
for and Fig. 8 and Fig. 9. 

The onset of the classical Coulomb phase is signalled by 
the sharpening of bow-tie like “pinch point” features around 
reciprocal lattice vectors such as q = (0,0,2). In the T ^ 0 
limit these features become sharp singularities^^’but at fl- 


0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 



Fig. 8. (Color line) Spin structure factor 4S(q) [Eq. (25)] in the antiferro¬ 
magnetic (Coulomb phase) region of the phase diagram [Fig. 2] of ‘Hbreathing 
[Eq. (1)], shown for ^ and = yg. With decreasing temperature, bow 
tie-like pinch point structures appear near certain Brillouin zone centers, e.g. 
q=(0, 0, 2). As r ^0 these become sharp singularities. At finite temperature 
the singularity is cut off by a correlation length which diverges as T~^l^, 
signalling the onset of algebraic (~ r“^) correlations in real space. The di¬ 
vergent part of the correlation length is proportional to the product JaH and 
therefore the development of the correlations is slower when closer to the 
isolated tetrahedron limit. 


nite tempeature they have a flnite width in the longitudinal 
direction The length scale ^\\ cuts off the algebraic real 
space correlations at long distances, reverting them to an ex¬ 
ponential form for r >> As T ^0, diverges as 
and the correlations are algebraic out to infinite separation. 

We may understand this crossover in the breathing py- 
rochlore case, by calculating as a function of Ja,Jb and 
T within the SCGA. This is done by expanding the SCGA 
prediction for <S(q) around (0,0,2). We And 

+ (43) 

The square of the correlation length may be written 
as a sum of two terms, one of which reflects single tetrahe- 






9 























J. Phys. Soc. Jpn. 


FULL PAPERS 



(a) 0 = fL; T= 5(»), I = 2{m), I = !(♦), I = 0.5(a), I = 0.2( ) 



(b) = If; £ = 5(»), £ = 2(.), £ = !(♦), £ = 0.5 (a), £ = 0.2( ) 


Fig. 9. (Color online) Angle integrated structure factor <Spow (q) [Eq. (27)] 
in the Coulombic spin liquid phase of ‘Hbreathing [Eq. (1)]. Correlations are 
shown for ^ and ^ = yf [cf. Fig. 2]. The development of sharp pinch 
point features, as seen in Fig. 8, is hidden by the angular integration and the 
powder structure factor shows a growing, broad peak near ^. This is consis¬ 
tent with experimental observations on the breathing pyrochlore LiInCr 408 
as discussed in Section 6 and shown in Fig. 10. 


0.6r 


0.5 



^‘8.0 0.5 1.0 1.5 2.0 

Q (A-b 


Fig. 10. (Color online) Comparison of SCGA calculations with experi¬ 
mental neutron scattering data^^ for a powder sample of for LiInCr 408 at 
T = SOAT. The lattice constant ao = 8.42A and the bond lengths dA = 
2.903A and ds = 3.052A are set to their experimental values.^^ The ex¬ 
change constants are taken from an estimate by Okamoto et al based on an 
empirical realtionship between the bond lengths and the exchange parameters 
in Cr spinels:^^ Ja = 60K and / 5 = 6 K. The SCGA calculation of the scat¬ 
tering has been multiplied by the form factor for Cr^+ ions \F{Q)\^. There 
is one adjustable parameter in the comparison which is a multiplicative scale 
factor setting the overall normalisation. The agreement is very good, which 
suggests that the estimated exchange parameters are approximately correct in 
this temperature range and underlines the validity of the SCGA description 
in the paramagnetic, cubic symmetry phase of LiInCr 408 . 


(48) 

^0 << ^Lulomb ^ ^11 ^ ^Coulomb ~ (49) 

and the correlation length diverges as as expected for a 
Coulomb phase. In the case Ja = Jb these expressions repro¬ 
duce the results of Conlon and Chalker.^^ 

In the high temperature limit we have 



dron correlations and dominates at high temperature the other 
of which reflects the correlations of the Coulomb phase and 
dominates at low temperature 

^f\ = ^0 + ^Coulomb (44) 

2 ^ ( ^pkt){Ja + Jb) \ 

^0 Uy V4d(r)2 + 8ySd(r)(/A + /fi)y ^ ^ 

^2 _ / 160^JaJb \ 

^Coulomb - \4A{Tf + m{T){JA + JB) J 

(46) 


In the spin liquid region of the phase diagram the Lagrange 
multiplier A{T) tends to a finite value at both high and low 
temperature. Thus, in the low temperature limit we have 


- (?)' 

y ^2 _ ( WJaJb 

r™ ^Coulomb - U y \A{T){Ja + Jb) 


(47) 


1 

T 


^Coulomb « ^0 ^ ^11 ^0 (52) 


and the correlations are essentially described by a model of 
isolated tetrahedra. Comparison of the length scales ^coulomb 
and ^0 gives a criterion to define the thermal crossover from 
isolated tetrahedron correlations to isolated Coulomb phase 
correlations in an antiferromagnetic breathing pyrochlore. 

In the powder structure factor shown in Fig. 9 the sharpness 
of the pinch point features is hidden by the angular integration 
[Eq. (27)] and the scattering simply shows a broad maximum 
around 


Qao 

2n 


= 2 . 


(53) 


This maximum grows in intensity and becomes sharper with 
decreasing temperature but remains always rather broad due 
to the angular integration. This is consistent with neutron scat¬ 
tering experiments on LiInCr 408 , as will be discussed further 
in Section 6. 
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6. Discussion 

In this article we have considered the minimal model for 
for magnetism on the breathing pyrochlore lattice [Eq. (1)]. 
We have seen that the ground state manifold depends only on 
the signs of the two exchange couplings, associated with the 
two species of tetrahedra [Section 2]. The classical ground 
state phase diagram [Fig. 2] contains a ferromagnetic phase 
(for Ja < 0, < 0) a disordered Coulomb phase (for Ja > 

0,Jb > 0) and a phase where an 0(L) degeneracy is lifted by 
the order-by-disorder mechanism (for Ja > 0, Jb < 0 or vice 
versa). The temperature development of the spin correlations 
in each case is calculated using the Self Consistent Gaussian 
Approximation (SCGA)"^^"^^ and is discussed for each region 
of the phase diagram in Sections 3- 5. 

At this point it is interesting to compare our calcula¬ 
tions with the results of the recent neutron scattering exper¬ 
iments by Okamoto et al?^ These authors measured the neu¬ 
tron scattering response of powder samples of breathing py- 
rochlores LiIn;cGai_;cCr408. Their neutron scattering data for 
X = \ (i.e. for LiInCr408) dXT = 30K, are shown in Fig. 10. 

In Ref. 30 the authors provided an estimate of the exchange 
parameters for FiInCr 408 , based on an empirical relationship 
between the exchange parameters and the bond distances, ob¬ 
taining 

Ja = 60K, Jb = 6K. (54) 

We have compared the results of our SCGA calculation us¬ 
ing this parameter set with the neutron scattering data from 
Ref. 32. The result is shown in Fig. 10. The comparison shows 
very good agreement, with one adjustable parameter which 
is a multiplicative scale factor setting the overall normalisa¬ 
tion. This suggests that the exchange parameters proposed in 
Ref. 30 [Eq. (54)] are at least approximately correct, at tem¬ 
peratures above the structural transition and underlines the va¬ 
lidity of the description presented here. If we fix the overall 

energy scale / = + /^ [Eq. (2)] and vary the ratio ^ 

we find that reasonably good agreement is obtained for the 
range of values 

0.05 < — < 0.15. (55) 

Ja 

Rather surprisingly, we are unable to obtain a convincing fit 
to the higher temperature data at T = 15OK. As noted by the 
authors of Ref. 32, the experimental scattering dXT = 150K 
is simply form-factor like, suggesting vanishing spin corre¬ 
lations, whereas our analysis predicts measurable nearest- 
neighbour spin correlations at T = 15OK for the estimated 
Ja = 60K and Jb = 6K. If we allow Ja and Jb to vary, we 
are unable to obtain any parameter set which simultaneously 
describes the data at both temperatures. 

This may be an indication that coupling to the lattice gives 
rise to an effectively temperature dependent exchange energy 
scale. This seems plausible since the dominant exchange path¬ 
way is direct exchange between Cr^+ ions, which is highly 
sensitive to the Cr-Cr distance. 

In conclusion, the recently synthesized breathing py¬ 
rochlore materials offer new opportunities in the study of frus¬ 
trated magnetism, playing host to spin-liquid and order-by¬ 
disorder physics. It may be particularly interesting for fu¬ 
ture experimental studies to investigate sulfide and selenide 


breathing pyrochlores^^ which could potentially realise the 
case of opposite exchange parameters discussed in Section 4, 
which exhibits an unusual structure of spin correlations and 
order-by-disorder. 

The determination of the full phase diagram of 
"^breathing [Eq. (1)] in the presence of quantum fluctua¬ 
tions remains an important and interesting open question. 
Moreover, experimental results on the breathing pyrochlores 
LiInCr 408 and LiGaCr 408 make it evident that the effect 
of coupling to the lattice is an important direction for 
future work.^^“^^ We expect that as in the “non-breathing” 
case,^^’^^’^^ the physics uncovered by a treatment of the 
spin-lattice coupling will be very rich indeed. 
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